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ABSTRACT 


Experimental  tests  have  been  made  with  various  configurations  of 
the  BALLUTE  in  order  to  develop  a  stable  retardation  device  for  meteoro¬ 
logical  rocxetsondes.  This  report  discusses  the  reduction  and  analysis  of 
these  tests.  Several  of  the  BALLUTE  configurations  are  shown  to  satisfy 
the  project  goal  of  providing  the  required  stability  as  well  as  a  sufficiently 
slow  fall  velocity  to  accurately  measure  winds  and  temperature. 
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SECTION  I 


INTRODUCTION 

Twc  types  of  BALLUTE  configurations  are  presently  being  tested 
and  evaluated  as  retardation  devices  for  Meteorological  Rocketsondes.  The 
two  configurations  are  the  8  or  12 -gore,  14-16  foot  diameter  BALLUTEs 
designed  for  use  with  the  Arcasonde  meteorological  package  and  the  light 
weight  4,  6,  or  12-gore,  6-7  foot  diameter  BALLUTEs  designed  co  be  used 
with  the  Loki-Dartsonde  package. 

The  mission  of  both  BALLUTE  configurations  is  to  provide  a  stable 
retardation  of  the  sonde  during  its  descent  from  rocket  apogee  (approximately 
200,  000  ft)  to  80,  000  feet.  Freedom  from  large  coning  angles  experienced 
with  the  more  conventional  silk  parachutes  will  result  in  cleaner  telemetry 
reception.  A  slower  fail  rate  is  also  desirable  since  it  will  allow  the 
BALLUTE  more  time  to  sense  the  atmosphere  and  thus  temperature  and 
winds  can  be  measured  to  a  higher  degree  of  accuracy.  In  order  to  attain 
the  .e  objectives  design  requirements  demand  that  the  BALLUTE  be  aero- 
dynamieally  a  very  stable  object  and  further  that  its  fall  velocity  be  sufficiently 
slow.  Sufficiently  slow  is  taken  to  mean  a  maximum  vertical  velocity  of 
300  ft/ second  at  180,  000  feet  or  a  ballistic  coefficient  value  of  W/C^A  =  0.  50. 
Slower  fall  rates  are,  of  course,  desirable. 


SECTION  II 


LOG  OF  FLIGHTS 

BALLUTE  Arcasonde  and  Dartsonde  flights  began  in  April  of  1966 
at  Cape  Kennedy.  A  log  of  the  succe-^ful  flights  and  their  characteristics 
are  given  in  Table  1, 
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BALLUTE  CHARACTERISTICS 
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SECTION  in 


DATA  REDUCTION 

The  data  from  the  successful  flights  were  reduced  at  UDRI.  The 
reduced  data  were  presented  in  the  form  of  a  blue  data  booklet  for  each  flight. 
The  format  of  the  booklet  and  what  it  contains  follow. 

Title  Page  --  Contains  launch  and  flight  log  information  as  well  as 
any  special  comments  concerning  performance  of  the  sonde  and  radars. 

Tabulations  --  As  many  as  five  (5)  sets  of  tabulations  may  be 
contained  in  a  booklet  if  all  temperature  information  is  available-  The  first 
tabulation  gives  the  component  velocity  and  acceleration  of  the  BALLUTE  as 
obtained  by  a  39-7  linear  linear  least  squares  fit.  The  second  tabulation 
computes  the  ballistic  coefficient  using  the  temperature  from  the  strip  chart, 
and  the  vertical  velocity  obtained  from  the  first  tabulation.  The  remaining 
tabulations  contain  the  rawintonde  flight  of  the  day  and  the  Cape  Kennedy 
BALLUTE  reduction  of  winds  and  density. 

Plots  —  If  temperature  data  is  available,  three  types  of  plots  are 
presented  in  the  booklets;  a)  vertical  velocity  v«>  aliitude,  b)  temperature 
vs  altitude,  and  c)  density  vs  altitude.  If  the  temperature  trace  is  poor  or 
not  available  only  the  vertical  velocity  vs  altitude  plot  is  included. 

Table  2  shows  the  content  of  each  data  booklet  that  was  reduced  by 

UDRI. 
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Cape  Kennedy 
Reduction 
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BALLUTE 


Flight 

Number 

Type 

of 

Flight 

Time 

Zulu 

Flight 

Date 

3084 

Ballute 

1606 

04-22-66 

2643 

Ballute 
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05-24-66 

2561 

Ballute 

1838 

05-26-66 

2793 

Ballute 

1445 

05-27-66 

0906 

Ballute 

1615 

07-06-66 

0715 

Ballute 

1509 

07-08-66 
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Parachute 
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07-08-66 
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08-15-66 

5554 
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1500 

08-16-66 
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04-12-67 
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05-09-67 
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05-10-67 

8622 
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05-10-67 

8767 

Ballute 

1527 

05-12-67 

6513 

Ballute 

1627 

05-12-67 

8534 

Ballute 

1115 

05-18-67 
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SECTION  IV 


METHOD  OF  REDUCTION 

4.  1  Velocity  and  Acceleration  of  Ballute 

The  velocity  of  the  BALLUTE  was  obtained  by  a  linear  least 
squares  fit  of  39-1/2  second  position  points.  The  slope  of  the  fitted  line  is 
called  the  velocity  at  the  midpoint  in  time  of  the  interval.  Tf  muar 
records  at  increments  less  than  1/2  second  then  thece  position  points  are 
averaged  to  obtain  1/2  second  increments  between  consecutive  coordinates. 

To  obtain  an  acceleration,  seven  (7)  velocities  are  fitted  linearly  and  the 
slope  assigned  to  the  midpoint  as  the  acceleration.  To  obtain  derivatives 
each  second  one  slides  two  position  coordinates  and  one  velocity  coordinate 
and  repeats  the  fitting  procedure.  For  a  further  description  see  Engler 
1965. 

4. 2  Temperature 

Temperature  was  either  read  directly  from  the  strip  chart  or  if  it 
had  already  been  reduced  it  was  interpolated  from  the  reduced  data  at  altitude 
increments  of  1000  feet. 

4.  3  Density 

The  density  was  calculated  from  temperature  by  using  the  Gas  Law 
and  Hydrostatic  Equation.  The  formula  is 

T 

<n(pn)  =  la  (pQ>  -  ^  J  j  dz  (See  Appendix)  (1) 

°  zo 


where  =  initial  density  -  obtained  from  rawinsonde  flight  for  the  day 

T  =  Temperature 

K  =  Gas  constant/molecular  weight  of  atmosphere 


=  gravitational  acceleration 


The  integration  started  with  the  lowest  BALLUTE  altitude  and  proceded 
upwards  so  that  the  initial  density  (pq)  could  be  taken  from  the  rawinsonde 
flight  for  that  day. 
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4.  4  Ballistic  Coefficient 

The  ballistic  coefficient  is  computed  from  the  equation 


where  a  =  ballistic  coefficient 
Z  =  vertical  velocity. 
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SECTION  V 


ANALYSIS 


5.  1  Wind 

Computation  of  horizontal  winds  has  not  been  included  in  the  present 
reduction  of  BALLUTE  flights.  The  following  analysis  is  designed  to  determine 
the  necessity  and  importance  of  the  acceleration  and  apparent  mass  of  the 
BALLUTE  in  a  wind  computation. 

Ignoring  the  Coriolis  force  and  the  buoyancy  of  the  payload,  the 
equations  of  motion  in  the  X  and  Z  direction  for  the  BALLUTE  can  be  written 
as: 


(m  +  m  )  X  =  i  p  C  AV  (W  -  X)  +  m*  (W  -  X) 
p  g  2  D  x  x 


(2) 


where 


(m  +  m  )  Z  =  I  pC_  AY  (W  -  Z)  +  m'  (W  -  Z)  -  m  |g|  (3) 

p  g  2  D  z  z  p 


irip  =  mass  of  BALLUTE  plus  payload 

m  =  mass  of  gas 
g 

m'  =  apparent  mass  of  BALLUTE 


C. 


Drag  Coefficient 


A  =  Cross  sectional  area  of  BALLUTE 

P  .  2  •  2*1  1^Z 

V  =  (X  -  +  (Z  -  w/ 


Assuming  the  absence  of  vertical  winds,  Equation  (2)  and  (3)  are  combined 
to  give  the  wind  equation: 


W-KW  =  X  -  K.  X 
x  2  x  1 


(4) 


where 


K.  = 


(m  +  m  +  m ')  Z 

“  Ji 


1  (m  +  m  +  m')  Z  +  m  |g|  ^2  (m  +m  +  m')  Z  +  m  |g| 


m '  Z 
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If  the  wind  equation  is  written  as  W  =  X  -  K,  X  +  K_  W  then 

x  1  2  x 

the  component  terms  on  the  right  hand  side  of  this  equation  can  be  plotted 
for  a  given  wind  field  to  determine  their  contribution  to  W  .  The  term: 

•  •  •  X 

X  represents  the  velocity  of  the  BALLUTE;  K^X  represents  the  contribution 
to  Wx  as  seen  by  the  acceleration  of  the  BALLUTE;  and  K2WX  represents 
the  contribution  to  Wx  due  to  the  apparent  mass  of  the  BALLUTE. 

Consider  a  Sinusoidal  Wind  Field 


Figures  1-4  are  plots  of  the  component  terms  of  the  wind  equation 

2irZ 

for  the  wind  field  W  =  A  sin  From  the  plots  it  is  seen  that  the 

X  fa  WV  V 

term  W^  contributes  less  than  1/2%  to  the  amplitude.  K^,  however, 
was  computed  using  m*  =  1/2  p  Vol  which  is  the  apparent  mass  formula  for 
a  sphere.  This  assumption  was  necessary  due  to  the  lack  of  knowledge 
concerning  the  apparent  mass  of  the  BALLUTE.  Even  if  the  formula  were 
considerably  in  error  the  apparent  mass  term  would  still  not  be  a  significant 
term  in  the  wind  equation.  The  acceleration  term  KjX,  however,  is 
significant,  particularly  at  high  altitudes.  It  is  also  seen  from  the  plots  that 
the  Dartsonde  BALLUTE  is  slightly  more  wind  sensitive  than  the  Arcasonde 
BALLUTE.  This  is  due  to  its  slower  fall  velocity  and  smaller  ballistic 
coefficient  (see  Table  3). 

Consider  a  Linear  Wind  Field 


The  sensing  ability  of  the  BALLUTE  can  also  be  analyzed  under  the 
influence  of  a  linear  wind  field.  Figures  5  and  6  compare  Arcasonde  BALLUTE 
velocity  to  wind  velocity  for  linear  winds  with  shear  values  from  .  013  to  .  1  per  sec. 
In  all  cases,  BALLUTE  velocity  is  not  an  accurate  measure  of  wind  velocity 
and  consequently  the  acceleration  of  the  BALLUTE  must  be  included  in  any 
wind  computation.  It  can  easily  be  shown  that  the  apparent  mass  term  is 
insign;ficant  for  a  linear  wind  field  as  well  as  for  a  sinusoidal  field. 

Hence  accurate  winds  can  be  determined  only  if  velocity  and 
accelerations  are  measured  accurately.  Initially,  velocity  and  acceleration 
have  been  obtained  from  the  39-7  linear  fit.  This  technique  was  used  as  a 
result  of  our  previous  experience  with  the  ROBIN  system  taking  into  con¬ 
sideration  the  slower  vertical  velocity  of  the  BALLUTE  by  lengthening  the 
position  smoothing  interval  from  31  to  39  points.  It  is  not  believed  that 
this  will  ultimately  be  the  smoothing  technique  decided  upon  to  compute 
velocities  and  accelerations  and  thus  calculate  a  wind.  'Indications  are 
that  the  vertical  smoothing  interval  should  not  necessarily  be  of  the  same 
length  as  the  horizontal.  Also  in  some  cases  a  cubic  velocity  fit  may 
improve  upon  a  linear  fit. 
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Figure  2.  Response  of  Dartsonde  BALLUTE  to 
Sinusoidal  Wind  Field 
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Figure  3.  Response  of  Arcasonde  BALLUTE  to 
Sinusoidal  Wind  Field 
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Figure  4.  Response  of  Dartsonde  BALLUTE  to 
Sinusoidal  Wind  Field 
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5.  2  Ballistic  Coefficient 


In  order  to  insure  an  acceptable  fall  rate  the  BALLUTE  has  a 
design  criterion  o:  maintaining  a  maximum  ballistic  coefficient  of 
W 

q  =  — — j  =  0.  050.  The  ballistic  coefficient  equation  is  derived  from  the 

cDA 

vertical  equation  of  motion  as  follows: 


(m  +  m  )Z  =  ip  C_  AV  (W  -Z)+m*(W  -  Z)  -  m  Igl 
p  g  2D  z  z  p  1  1 

If  over  a  sufficiently  small  interval  Z  can  be  considered  constant  and 

W  =0  the  equation  simplifies  to: 
z 


"V  1*1 

cda 


'2  P 


V  z. 


(5) 


The  common  formula  for  ballistic  coefficient 


q 


“pi*1 

V 


(6) 


is  obtained  by  further  assuming  that  V  =  Z.  This  assumption  is  valid 
except  for  the  first  3000  feet  of  descent.  During  this  period  the  BALLUTE 
has  a  large  horizontal  velocity  transmitted  to  it  from.the  rocket,  and  thus  its 
velocity  relative  to  the  air  may  be  much  larger  than  Z. 


In  the  data  booklets  two  methods  have  been  used  to  evaluate  q. 


Method  1 


From  the  thermistor  temperature  trace  the  density  is  computed 
by  F.quation  (1),  Z  is  obtained  from  the  39-7  linear  fit,  and  q  *s  then 
1  •  2 

computed  as  q  =  -  p  Z  .  A  ballistic  coefficient  is  computed  by  this  method 
at  each  1000  feet  of  altitude. 

Method  2 


The  second  method  of  calculating  a  ballistic  coefficient  assumes 
that  the  BALLUTE  fell  in  the  1962  Standard  Atmosphere.  For  a  given 
ballistic  coefficient  (q )  a  Z  vs  Z  curve  is  generated  from  Equation  (3)  by 
obtaining  density  as  a  function  of  altitude  from  the  1962  Standard  Atmosphere. 
A  series  of  these  Z  vs  Z  curves  are  generated  for  values  of  q  from  q  =  0.  01 
to  q  =  0.  10. 
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From  a  reduced  BALLUTE  flight  a  plot  is  made  of  Z  vs  Z.  This 
velocity  profile  of  the  BALLUTE  is  then  compared  to  the  ballistic  coefficient 
curves.  The  curve  that  most  closely  follows  the  BALLUTE  velocity  profile 
is  considered  a  representative  value  for  q.  It  has  been  experimentally 
varified  by  both  Method  1  and  Method  2  that  the  value  of  q  is  relatively 
constant  over  the  entire  flight.  If  the  atmosphere  is  not  model  an  error  will 
be  made  using  Method  2  to  compute  a  ballistic  coefficient.  For  example  if 
density  is  only  80%  of  Model,  then 


1 


.  2 

q  =  -  (0.  8  p  )  Z 

"True  2  1  PModel' 


.  2 

q  =  -  (P  )  Z 

^Method  2  2  'KModel 


1 


% 


^True  ^Method  2 


Error 


‘True 


?  <°-  2  ^Model>  ^ 
?  (°-  8  PModel>  ^ 


=  25% 


Table  3  gives  the  ballistic  coefficient  for  each  flight  as  determined  by  both 
methods.  The  12-gore  Dartsonde  BALLUTEs  show  a  somewhat  smaller 
ballistic  coefficient  (.  045)  than  the  12-gore  Arcasonde  BALLUTEs  (.  05). 
However,  a  later  series  of  firing  using  a  4-gore  Dartsonde  BALLUTE  and 
8-gore  Arcasonde  BALLUTEs  have  produced  a  substantial  decrease  in  the 
ballistic  coefficient.  An  average  value  of  ballistic  coefficient  for  the  4-gore 
Dartsonde  BALLUTE  is  0.  035  and  the  8-gore  Arcasonde  BALLUTE  0.  04. 

These  ballistic  coefficients  provide  the  required  descent  rate. 

5.  3  Temperature 

In  order  to  insure  an  accurate  measure  of  temperature,  a)  the 
vertical  velocity  of  the  BALLUTE  must  be  sufficiently  slow  to  insure  that 
induced  aerodynamic  heating  is  negligible  and  b)  the  BALLUTE  is  aerody- 
namically  a  stable  vehicle  so  that  the  system  does  not  oscillate  or  come 
through  wide  angles  thus  resulting  in  aerodynamic  heating  of  the  thermistor 
as  well  as  signal  dropouts.  The  established  vertical  velocity  criterion  is  to 
maintain  a  ballistic  coefficient  <  0.  05.  As  can  be  seen  by  Table  3,  BALLUTEs 
achieve  this  result  in  nearly  all  canes.  The  stability  condition  has  been 
verified  by  an  observation  of  position  plots  as  well  as  the  temperature  strip 
chart  which  shows  very  few  signal  dropouts.  (See  Wright  and  Graham,  1966.  ) 
Figures  7-12  are  North-East  position  plots  of  flight  906.  Throughout  the 
entire  flight  path  there  is  no  evidence  of  any  periodic  c.erodynamic  oscillation 
in  the  BALLUTE's  path.  The  very  smooth  path  prevalent  over  several  large 
segments  also  substantiates  the  stability  of  the  BALLUTE. 
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TABLE  3 


BALLISTIC  COEFFICIENTS 

Flight  Ty  pe 

Number  of  Flight 

Time 

Zulu 

Flight 

Date 

Ballute 

Series 

Ballistic 

Coefficient 

Method  1 

Method  2 

3084 

Ballute 

1606 

04-22-66 

Dart  2 

.  05 

2643 

Ballute 

1500 

05-25-66 

Dart  3 

.  035-. 

045 

.  045 

2561 

Ballute 

1838 

05-26-66 

Dart  4 

.04  -. 

05 

.  045 

2793 

Ballute 

1445 

05-27-66 

Dart  5 

.06  -. 

07 

-065 

0906 

Ballute 

1615 

07-06-66 

Dart  7 

.  035-. 

04 

.  04 

1715 

Ballute 

1509 

07-08-66 

Areas  6 

.04  -. 

07 

.  04 -.  06 

7683 

Parachute 

1630 

07-08-66 

.05-.  07 

7391 

Ballute 

1500 

08-15-66 

Areas  9 

.08  -. 

09 

.  08 

5554 

Ballute 

150C 

08-16-66 

Areas  10 

.055-. 

065 

.  06 

6097 

Parachute 

1600 

08-16-66 

.04  -. 

05 

.  055 

7867 

Ballute 

1500 

08-18-66 

Areas  11 

.05  -. 

06 

.055 

2419 

Parachute 

1359 

09-14-66 

.  04 

2691 

Ballute 

1830 

09-14-66 

Areas  12 

.05  -. 

06 

.05 

1615 

Parachute 

1359 

09-15-66 

.05 

3163 

Ballute 

1600 

09-15-66 

Areas  13 

.05 

.05 

1947 

Ballute 

1500 

09-19-66 

Areas  14 

.045-. 

05 

.05 

1731 

Ballute 

1500 

09-20-66 

Areas  15 

.05 

.05 

3309 

Parachute 

.05 

4765 

Ballute 

1900 

11-14-66 

Dart  11 

.04 

3166 

Ballute 

2007 

11-14-66 

Dart  12 

.  045 

6258 

Ballute 

1520 

11-15-66 

Dart  13 

.035 

4287 

Ballute 

1920 

11-17-66 

Areas 

.05  -. 

06 

.05 

extra 

3729 

Ballute 

1830 

11-18-66 

Areas 

.05 

extra 

2577 

Ballute 

1430 

01-30-67 

Areas  16 

.  06 

0015 

Ballute 

1703 

01-30-67 

Dart  15 

.05 

2694 

Ballute 

1500 

01-31-67 

Areas  17 

.  055 

0099 

Ballute 

1608 

01-31-67 

Dart  16 

.  035 

1616 

Ballute 

1500 

02-27-67 

Dart  17 

.  035 

0613 

Ballute 

1630 

02-28-67 

Dart  18 

.035 

0776 

Ballute 

1500 

03-01-67 

Areas  18 

.015 

7092 

Ballute 

1600 

03-02-67 

Areas  19 

.01-.  1 

3605 

Ballute 

1500 

04-12-67 

Areas  20 

.04  -. 

05 

.  04-.  055 

0739 

Ballute 

1500 

04-19-67 

Areas  21 

.04  -. 

05 

.05 

5395 

Ballute 

1421 

05-09-67 

Dart  21 

8990 

Ballute 

1500 

05-10-67 

Areas  22 

.04 

.045 

8622 

Ballute 

1600 

05-10-67 

Areas  23 

8767 

Ballute 

1527 

05-12-67 

Dart  23 

.  025-. 

03 

.  025 

6513 

Ballute 

1627 

05-12-67 

Dart  24 

.  025-. 

03 

.  03 

8534 

Ballute 

1115 

05-18-67 

Dart  25 

.  035 

.  035 
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Figure  7.  North  vs  East  Position 


Figure  8.  North  vs  East  Position 


Time  increasing  from  right  to  left 
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Figure  11.  North  vs  East  Positron 


Figure  12.  North  East  Position 


5. 4  Density 


Density  is  obtained  from  temperature  by  use  of  Equation  ().).  The 
integral  is  evaluated  by  use  of  the  trapezoidal  rule  (see  Appendix).  Integration 
by  the  trapezoidal  rule  assumes  a  linear  interpolation  of  temperature  between 
data  points,  which  for  the  purposes  of  reduction  were  read  every  1000  feet. 

The  initial  density  required  for  Equation  *1)  is  obtained  from  the 
rawinsonde  flight  for  the  day  of  launch.  If  the  initial  density  is  in  error  by 
eg  this  will  introduce  a  constant  percentage  error  in  density  for  the  entire 
flight  given  by: 

e 

%  Error  =  —  .  (See  Appendix.  ) 
p0 


SECTION  VI 


CONTINUING  INVESTIGATION  OF  BALLUTE  AS  WIND  SENSOR 


Future  investigation  at  UDRI  will  be  directed  towards  determining 
an  optimum  method  to  compute  winds  and  the  resulting  error  in  winds  as 
a  function  of  radar  error  using  this  method.  The  UDRI  is  presently 
developing  and  evaluating  smoothing  methods  for  obtaining  velocity  and 
acceleration  using  the  criterion  of  a  minimization  of  the  total  error  in  the 
wind.  The  total  error  is  defined  as:  <r2.  =  <r„,,  +  (r£,  .  where  cA.  , 


Fit  Noise 


Noise 


is  the  random  error  in  wind  due  to  the  noise  in  the  position  coordinates  and 
(rj'.j  is  the  bias  error  in  wind  introduced  by  th*  lack  of  fit  of  the  smoothing 


function.  The  UDRI  also  is  investigating  the  effect  of  the  radar  range 
resolver  error  on  wind  determination.  Considerable  work  has  already 
been  done  with  regards  to  these  two  aspects  on  the  High  Altitude  ROBIN, 
and  ROSE  systems.  Much  of  this  work  can  be  applied  to  the  BALLUTE 
system. 
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APPENDIX 


DENSITY  OBTAINED  FROM  A  TEMPERATURE  PROFILE 


I. 


1)  p  -  K  p  T  Gas  Law 

2)  dp=  -  p  jg  |  dZ  Kydrostauc  Equation 

where 


p  =  pressure 
p  =  density 
T  =  Temperature 
g  =  gravitational  acceleration 
Z  =  Altitude 

K  =  Universal  Gas  Constant/Molecular  weight  of  atmosphere 
From  Equation  (1) 

dp  =  K  p  dT  +  KT  dp 
substituting  into  Equation  2  we  get 

-  p  jg  |  dZ  =  K  p  dT  +  KT  dp. 

Simplifing 


Ip  -  isL 

p  KT 


dZ 


dT 

T 


Initial  Conditions: 


A»  z  =  z0,  g  =  g0,  T  =  T0.  p  =  p0 

AtZ  =  z,g  =  g,T  =  T,p  =  p 

n  n  n  n 

Integrating  down  from  the  initial  altitude  Z^we  get; 
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or 


T 

Ln  (p  )  =  Ln  (p  )  -  Ln  (-H) 
n  u  i. 

0 


1M 

T 


dZ 


Assuming  a  linear  interpolation  of  temperature  between  data  points,  the 
integral  is  evaluated  by  the  trapezoidal  rule. 


Ln  (pn)  =  Ln  (pQ)  -  Ln  )  - 

0 


1 

2K 


gn  gl  gn  1  gn 

<rn  +  r"zi-zo>  +  -  M^^xzn-zn 

0  1  n-1  n 


II.  Error  Due  to  the  Initial  Guess  of  p 


Let  eg  be  the  error  in  the  initial  guess  of  pg.  Then  the  computed 
value  of  pn  is  obtained  from; 


Ln  (p  )  =  Ln  (p  +  e  )  -  Ln  { 
n  u  u 


p„  =  <P0  +  V  (f“  1  EXP  ( 

n 


To  KJZ 

0 

fz" 

.1  IsL 

KJz.  T 


U1 


dZ 


dZ) 


The  percentage  error  is  given  by: 


%  Error 


|PnTrue  Pn  Comp| 
Pn  True 


e0(^)E*P 


n 


<■£ 


M 

K 


dZ) 


t  .  r  z 

f/)Exp< -R  “ 

n  Jzo 


|l  dZ) 


%  Error 


a  constant  for  the  entire  flight. 
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